Many-Body Theory of Dilute Bose-Einstein Condensates with Internal Degrees of 

Freedom 
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The Bogoliubov theory of weakly interacting bosons is gen- 
eralized to Bose-Einstein condensates with internal degrees of 
freedom so that a single effective Hamiltonian produces vari- 
ous many-body ground states or metastable spin domains and 
the corresponding collective modes on an equal footing. 

PACS numbers: 12.20.Ds, 42.50.Ct, 42.50.Lc 



The internal degrees of freedom of alkali Bose-Einstein 
condensates (BECs) arise primarily from electronic spin 
and are hence by far more amenable to manipulation 
than those of superfluid ^He. The realization of "spinor" 
BECs has motivated many researchers to study the 
magnetism of superfluid vapors The excitation 

spectra of spinor BECs have been examined by Ohmi 
and Machida (2) and by Ho ||] using mean-field theory. 
The many-body ground states were investigated by Law 
et al. [Q and by Koashi and Ueda ||] and Ho and Yip ||] 
in the absence and presence, respectively, of a magnetic 
field; in the latter case, it is predicted that the magnetic 
sublevel m = of a spin-1 antiferromagnetic BEC be- 
comes populated as the magnetic field decreases. The 
main aim of this Letter is to present a unified theory 
that derives the various many-body states of BEC and 
the corresponding excitation spectra on an equal foot- 
ing. Previous work has been unable to address this 
problem, as it has been assumed that the internal degrees 
of freedom are independent of other degrees of freedom 
such as the density and spin waves. 

We consider a system of N spin-1 identical bosons in- 
teracting via a binary, contact interaction, subject to a 
uniform magnetic field B. The Hamiltonian of the sys- 
tem is given by 



H ^ dr 



-^V^+C/(r)) 



(1) 



where p{> 0) is the product of the gyromagnetic ra- 
tio and the external magnetic field B applied in the z- 
direction; f^^ = {fa^Jap^ faff) represents the compo- 
nents (a,/3 = 1,0,-1) of spin-1 matrices; and ^Q,(r) de- 
notes the field operator that annihilates a boson in mag- 
netic sublevel a at position r. The quadratic Zeeman ef- 
fect will be taken into account when it becomes relevant. 
The coupling constants cq and ci are related to the spin- 
singlet and triplet s-wave scattering lengths Os and at by 



Co — AttH^ {as -\- 2at) / 3M and ci — 47r?i^(at — as)/3M, re- 
spectively. In Eq. (|l|) the repeated indicies are assumed 
to be summed. 

When the system has translation invariance, it is 
convenient to expand the field operator as ^'^(r) = 
(l/\/y) flkae*'''", where V is the volume of the sys- 
tem. Equation (Q) then reduces to 



H 



-^Yl :(coPkPk + CiS[-Sk):, 



it 



(2) 



h k^/2M, nka = flk^aka, Pk 
k = LqQ/j fa/3«qa«q+k/3, and :: denotes 



where ek 

Sq^qa^q+ka, S 

normal ordering. 

When the interparticle interactions are repulsive, BEC 
occurs in the k = state. We therefore separate the 
k = components aoa in the Hamiltonian and keep 
the terms only up to second order in aka (k 7^ 0). Be- 
cause our interest is to derive various many-body states 
and the corresponding low-lying excitation spectra on an 
equal footing, wc should not replace the operators ooa 
by c- numbers. After some algebraic manipulation, we 
obtain H 

^eff ^ 'J^NiN - 1) - |liti - Pino, no^i) 



J2 J2 [ek-pa 

k/OQ=0,±l 



V 



-no- 



k^O 



+ Ci[ajj.^a^j^_^aoi-aaoa-i + akfj_iaLkQ'^OQ-iaoQ 

+ (^kQ"ka-i + '2al^_^ak^a)ao^-iaoa] + h.c.j (3) 
where h.c. denotes the hermitian conjugates of all 
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preceding terms in the curly brackets, = (a, 
2aQj^aQ_j)/V3 creates a pair of spin-singlet bosons when 
acting upon the vacuum, and -Dk = X]q ^oa'^ka and 
S'k = J2a Q^floa^ka dcscribc density-wave and spin-wave 
operators, respectively. The first three terms on the 
right-hand side of Eq. were studied in Refs. [^|-^; 
the remaining terms show how the internal and spatial 
degrees of freedom couple to one another. 

To set a reference frame for later discussions, let us 
recapitulate the main results of Refs. |^~^. The Hamil- 
tonian discussed in these references reads essentially 



1 



p{noi - no-i). 



(4) 



When ci < 0, the energy minimum is obtained by putting 
ah particles in the state with k = and a = 1: 



\N2 ^0,F = N,F^ 



N) = 1^2g^|vac;, 



(5) 



where F and F^ denote the total spin and its projection 
on the z axis, and N2 is given in terms of F and the total 
number of atoms N hy N2 = {N—F) /2. Because all spins 
are aligned in the same direction, the system is ferromag- 
netic. When ci > 0, the first term on the right-hand 
side of Eq. (^) energetically favors the spin-singlet (or 
antiferromagnetic) correlation, whereas the second one 
favors the parallel-spin configurations. For a given exter- 
nal magnetic field (oc p) , the energy minimum is attained 
for the state [|| 



\N2 - {N-F)/2,F,F,=F)„ 



vac) 



(6) 



where F denotes an integer that satisfies — 1/2 < 
p/ci <F + 3/2, and F_ = (/^ - ify)c.pal^aop. In both 
Eqs. (^) and the depletion of the condensate due to 
interactions is not taken into account. 

When the components with k = are macroscopi- 
cally occupied, we may replace aoa in Eq. (^) with the 
c- numbers a/TVCo- To find the correct excitation spectra, 
it is crucial to take the depletion of the condensate into 
account. This is done by setting 



rika 



(7) 



where the last term describes the fraction of the deple- 
tion. The energy of the condensate is given by 

Eo = -^iCo' - 2C1C-1P -pA^dCiP - IC-iP), (8) 

where (^a are determined by requiring that Eq be mini- 
mized subject to constraint (|^). 

For ci < 0, the minimum of Eq is reached when 

Co = c-i = o, ici|2 = i-l J2 ^ko. (9) 



The Hamiltonian (Bh then reduces to 



2V 

Ek^kl 



{cQ + ci)N 
2V 



+ (e + p)nko + {e + 2p —)nk-i 



(flkia-ki + Hki + h.c.) 



(10) 



Diagonalizing this Hamiltonian, we find the excitation 
spectra as 



i^k,i Vek[ek + 2(co + Ci)n], 
^k,o = ek + P, 
^l-i = ek + 2p- 2cin, 



(11) 



in agreement with Ref. j2j. We may gain some insight 
into the nature of the many-body correlation by writing 
the many-body wave function in the coordinate represen- 
tation: 

^'(riai,---,rjvaAr) = (vac|*ai (ri) (ri)|$). (12) 

Here 1$) is given in the Bogoliubov approximation by 

1$) - exp (/)oa|)i - ^ i^/ca[ia^ki ) l^ac), (13) 

where 0q = iV[l - (8/3)^na^/7r] and z/^ = 1 + Ck - 
\/ Ckick + 2) with Ck = k'^ / {STTatji) . In the dilute limit, 
the many-body wave function becomes very small unless 
all Qffe's are equal 1, and we find that ||] 

exp(-|(37r - 8)iV^naf/7r) 
~ (27riVF2^)i/4 



(14) 



where = jr^ — r^l and the terms that are of the order 
of 1/A^ are ignored. This result clearly shows that two 
bosons strongly repel each other when their distance be- 
comes smaller than the healing length ^ = {8'Katn)~^^'^ . 

For ci > 0, minimizing F = Eq — SJ2a ICap! where S 
is a Lagrange multiplier, yields 




ICol e 



i(</)i+0_i-20o) 



(15) 



where = argC^ and 7 = 2p/(cin). In the presence of 
the external magnetic field (i.e., 7^0), the minimum of 
F is attained when 



Co=0, IC±lP = ^ ± I - ^ E 



(16) 



and the corresponding Hamiltonian reads 

^ (ck + cinSa,Q)nwa + ^ E {<^o(^k-D-k + -D^^^k) 
+ci{SkS-k + SlSk + ?^|ClC-l|ak,oa-k,o) + h.c.j . (17) 
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The a = mode is decoupled and its dispersion relation 
is given by 



(18) 



The a = ±1 modes are coupled, and the dispersion rela- 
tions of the coupled modes are given by 



e\+ 2(co + ci)nek± ekW rt^(co-ci)2+ 



Cl 



again in agreement with Ref. 1^ . Our theory thus repro- 
duces all known collective modes that are derived using 
a different method. 

An interesting situation arises in the absence of a mag- 
netic field (i.e., 7 = 0), where condition ([T^ ) can be met 
if 



(20) 



This result is in accordance with a rule pointed out in 
Ref. |9j, that is, when the interaction is attractive, the 
relative phase coherence, as implied by the constraint 
(pO|), will spontaneously emerge if more than two BECs 
coexist. In the present case, the attractive force applies 
among three spin components, as implied by the condi- 
tion Cl > (see Eq. (^). 

When the relation ( pO| ) holds, Eq. ( p^ ) imphes that 
ICol is arbitrary, so we have |Ci| = |C-i|, ICoP = 
1 - 2|Cip. This, combined with Eq. (|^), leads to 
a vectorial order parameter as C = (Ci,CojC-i) = 
gi0o(_ei('^o-0-i) sin/3/V2, cos/3, e"*('^«""*-i) sin/3/V^). 
This result agrees with Eq. (5) of Ref. ^ and implies that 
the spin components can change under spatial rotation. 
It is surprising that the many-body ground state predicts 
otherwise: from Eq. (^, we find that at zero magnetic 
field (hence F — Q) spin components cannot change un- 
der spatial rotation but must always be the same [^|J^] 
and is given by fiQa = {N ~ ^^ka ^ka) /3 for a = 0, ±1, 
where we take into account the depletion of the conden- 
sate X^ka ""ka duc to interactions. This is a consequence 
of the fact that the true ground state is composed of spin- 
singlet "pairs" and is therefore invariant under rotation. 
The above result for rioa and the relative phase relation 
( 2(j) may be used to drastically simplify the Hamiltonian 



givmg 



rrAF 

-"7=0 



+ ^ [(ek+ Con)dl.d]^ + {e^ + cin)(sj^Sk + q^qk) 

k/O 



2V ^ ' 2V 



+ 2(corfkrf-k + ciSkS_k + ci^kg-k + h-c.) 



(21) 



dk = ^(Skie 
\/3 

Sk = -^(flkie"*"^' 
9k = -^(akie^^'^i 



akoe ^" 4 
flk-ie ^ 
2akoe-'^'' 



flk-ie 



flk-ie 



(22) 



We thus find that the density, spin, and quadrupolar fiuc- 
tuations provide independent excitations. The novelty of 
these modes is that they are phase-locked to the spin 
components of the condensate. Diagonalizing the Hamil- 
(19) tonian (^), we obtain the following excitation spectra: 



E^l = Vek[ek + 2co?i], 



77>AF 



pAF 
^q,k 



Vek[ek + 2cin]. 



(23) 



In the thermodynamic limit, these collective modes 
survive only at exactly zero magnetic field, and may 
be viewed as singular. However, in mesoscopic situa- 
tions, e.g., when the system is confined in a quasi-one- 
dimensional torus (to which the present theory applies 
with minor modifications), these modes survive at small 
magnetic fields, provided the k = 0, a = mode is 
macroscopically occupied (for a more precise definition 
of "small" magnetic fields, see Ref. §). 

Recently, an MIT group has observed the formation 
of metastable spin domains [0 . They first prepared all 
atoms in the a = 1 state and then placed half of them 
in the a = state by irradiating the rf field. Letting the 
system evolve freely while using the quadratic Zeeman 
effect to prevent the a = — 1 component from appearing, 
they found that spin domains formed with the two com- 
ponents alternatively aligned. It was discussed llO|-|l^] 
that this phenomenon is due to the imaginary frequen- 
cies of the excitation modes. Here we use our theory to 
confirm this hypothesis and to derive a general expression 
for the dispersion relation. The experimental conditions 
in effect amount to setting 



ICol 



^5](riki+riko), C-i=0. (24) 



Our Hamiltonian (W reduces to 



H = 'Jl±^niN - 1) - '-^nN 



E 

k^O 



Cqh 3cin\ / cqu cifi^ 
Ck + — + I nki + (^Ek + ^ ^ ' "ki 



Co 



+y'T-(akoa-ko + h-c.) 

Co + Cl ^ ^ ^ . ,i 

H ^ '^(akia-ki + flkia-ko + aki"-kO + h.c. 



(25) 



where 



Diagonalizing the Hamiltonian, we find the dispersion 
relations as 
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(Ek)^ =el + {2u + i;)ek + ± [4(w^ + 2uv + 2v'^)el 



(26) 



+2w'^(2u + v)ek -V^ [u + -V 



where u = CQn/2 and v = cin/2. For parameters of 
the MIT experiment n ^ 10^^ cm^^, at 29 A, and 
as ~ 26 A we find that w/w < 1. Ignoring in Eq. (|2^) 
higher-order powers of v/u, we obtain 



(El, 



el 



{2u + v)ek±2{u + v)ek, 



where the plus sign corresponds to the density wave, 
while the minus sign corresponds to the spin wave. We 
see that the energy of the spin wave becomes pure imagi- 
nary for Ek < V, implying the formation of spin domains. 
The corresponding wavelength defines the characteristic 
length scale of the spin domains as 



{at - as)n 



(27) 



This result agrees with that of Ref. except for a nu- 
merical factor. Using the above parameters, we obtain 
A ~ ISfim, in reasonable agreement with the observed 
value of about 40/im [l0|] . 

In conclusion, we have presented a versatile many-body 
theory of spin-1 dilute bose gas. Depending on the pa- 
rameters, a single effective Hamiltonian (||) describes var- 
ious many-body ground states, metastable spin domains, 
and the corresponding excitation spectra. An extension 
to spin-2 BEC can be similarly carried out and will be 
reported elsewhere. 
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